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ABSTRACT 

Weak gravitational lensing has become a common tool to constrain the cosmological 
model. The majority of the methods to derive constraints on cosmological parameters 
use second-order statistics of the cosmic shear. Despite their success, second-order 
statistics are not optimal and degeneracies between some parameters remain. Tighter 
constraints can be obtained if second-order statistics are combined with a statistic that 
is efficient to capture non-Gaussianity. In this paper, we search for such a statistical 
tool and we show that there is additional information to be extracted from statisti- 
cal analysis of the convergence maps beyond what can be obtained from statistical 
analysis of the shear field. For this purpose, we have carried out a large number of 
cosmological simulations along the a^-Vtrn degeneracy, and we have considered three 
different statistics commonly used for non-Gaussian features characterization: skew- 
ness, kurtosis and peak count. To be able to investigate non-Gaussianity directly in the 
shear field we have used the aperture mass definition of these three statistics for dif- 
ferent scales. Then, the results have been compared with the results obtained with the 
same statistics estimated in the convergence maps at the same scales. First, we show 
that shear statistics give similar constraints to those given by convergence statistics, 
if the same scale is considered. In addition, we find that the peak count statistic is the 
best to capture non-Gaussianities in the weak lensing field and to break the ag-^m 
degeneracy. We show that this statistical analysis should be conducted in the con- 
vergence maps: ffist, because there exist fast algorithms to compute the convergence 
map for different scales, and secondly because it offers the opportunity to denoise the 
reconstructed convergence map, which improves non-Gaussian features extraction. 
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1 INTRODUCTION 

Gravitational light deflection, caused by large scale structure 
along the line- of- sight, produces an observable pattern of 
alignments in the images of distant galaxies. This distortion 
of the images of distant galaxies by gravitational lensing, 
called cosmic shear, offers an opportunity to directly probe 
the total matter distribution of the Universe, and not just 
the luminous matter. Therefore, the statistical properties of 
this gravitational shear field are directly linked to the statis- 
tical properties of the total matter distribution and can thus 
be directly compared to theoretical models of structure for- 
mation. Despite some systematics (PSF distortion, intrinsic 
alignments...), this approach is extremely attractive since it 
is unaffected by the biases characteristic of methods based 
only on the light distribution. 
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Since its first detection (Van Waerbeke et al.l 2000 



'Kaise r et al.|2000|[Wittman etal.|2000| [Bacon et al.|200"of , 
cosmic shear has rapidly become a major tool to constrain 
the cosmological model (for review, see e.g, 



Mellier 



1999 



Bartelmann and Schneider|[200T] |Refregier||2003| |Hoekstra| 



and Jain||2008| [Munshi et al.|2008 ) 

In most weak lensing studies, second-order statistics are 
the most commonly used statistical probe 



(e.g. 



Maoli et al. 



2001] [Hoekstra et al.||2006l [Benjamin et al.|[20'07| |Fu et al 



2008 ) because of their potential to constrain the power spec- 



trum of density fluctuations in the late Universe. However, 
second-order statistics are not optimal to constrain cosmo- 
logical parameters. For example, they only depend on a de- 
generate combination of the amplitude of matter fluctua- 
tions (78 and the matter density parameter Q.rn { Maoli et al 



2001 Refregier et al.„2002, Bacon et al.„2003, Massey et al. 



20051 |Dahle[[2006t . 

The optimality of second-order statistics to constrain 
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cosmological parameters depends heavily on the assump- 
tion of Gaussianity of the field. However the weak lensing 
field is composed, at small scales, of non-Gaussian features 
such as clusters of galaxies. These non-Gaussian signatures, 
which can be measured via higher-order moments, carry ad- 
ditional information that cannot be extracted with second- 
order statistics. Since the non- Gaussianity is induced by the 
growth of structures, it holds important cosmological infor- 
mation. Many s tudies (e.g.^Bernardeau et al._1997a^^Takada 
and Jain|[2004l [Kilbinger and Schneider||2o65| |Pires et ~ 



2009a! 'Berge et al.|2010| have shown that combining second- 
order statistics with higher-order statistics tighten the con- 
straints on cosmological parameters. 

Most non-Gaussian studies (e.g. Schneider et al.||1998 
[KTlbinger and Schneider||2005| |Dietrichl 



Jar vis et al 



and Hartlap 



2004 



2010 



have been performed in the shear field 
because it can be directly derived from the shape of galax- 
ies. This paper aims to produce evidence that there is addi- 
tional information to be extracted from a higher-order sta- 
tistical analysis of the convergence maps beyond what can 
be obtained from a higher-order statistical analysis of the 
shear field because higher-order statistics are probing the 
non-Gaussian features of the signal and these non-Gaussian 
structures can be better reconstructed in convergence maps 
using a denoising. In Pires et al. (2009a), the efficiency of 



several higher-order convergence statistics have been com- 
pared to discriminate cosmological models along the ag-^m 
degeneracy. In this paper, we are interested in showing the 
advantage of using these higher-order convergence statistics 
compared to higher-order shear statistics. This comparison 
cannot be performed directly because the evaluation of non- 
Gaussian statistics in the shear field requires to use their 
aperture mass definition fSchneider et al.|[l998 ) or another 
different filter that is defined for a given scale 0. A fair com- 
parison with convergence statistics requires the statistics in 
the convergence maps to be estimated at the same scale. A 
stationary wavelet transform, the "a trous" wavelet trans- 
form, has been used in this paper to compute the conver- 
gence statistics at the given scale 0. 

The paper is organized as follows. In section 2, the cos- 
mological models selected for this study are described, fol- 
lowed by a short description of the weak lensing simulations. 
Section 3 summarizes the different statistics used in this 
study. We give the definition of the aperture mass Map and 
present the three shear statistics considered in this study. 
Then, the "a trous" wavelet transform is defined, as well as 
the three related convergence statistics. Section 4 presents 
our results and we summarize our conclusions in section 5. 



2 SIMULATIONS OF WEAK LENSING MASS 
MAPS 

N-body simulations have been used to numerically compute 
the variation of the different statistics with cosmologi- 
cal parameters and then compare their ability to break 
degeneracies. We carried out N-body simulations for 5 
different cosmological models along the ag-^m degeneracy 
corresponding to current constraints from a power spectrum 
analysis. The N-body simulations were carried out using 
the RAMSES code ( |Teyssier| [2002t . Fig. [l] shows the 
distribution of these cosmological simulations in the ag-^m 
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Figure 1. Location of the 5 simulated cosmological models in 
the (jg-^m plane. The 5 cosmological models have been selected 
along the crg-^m degeneracy corresponding to current constraints 
from a power spectrum analysis. 



plane. The characteristics of these cosmological models have 
been given in Pires et al. (J2009a). For each cosmological 
model, we have run 100 realizations in order to quantify 
the observational uncertainties. 



In the N-body simulations that are commonly used in 
cosmology, the dark matter distribution is represented by 
discrete massive particles. The simplest way of treating these 
particles is to map their positions onto a pixelized grid. In 
the case of multiple sheet weak lensing, we do this by taking 
slices through the 3D simulations. These slices are then pro- 
jected into 2D mass sheets. The effective convergence can 
subsequently be calculated by stacking a set of these 2D 
mass sheets along the line of sight, using the lensing effi- 
ciency function. This is a procedure that was used before by 



Vale and White ( 2003 ) , where the effective 2D mass distri- 



bution Ke is calculated by integrating the density fluctuation 
along the line of sight. Using the Born approximation, which 
assumes that the light rays follow straight lines, the conver- 
gence can be numerically expressed by 



SHQQrnL -^^ Xi{Xo 



Xi 



2c2 



Xoa{xi 



Nts^ 



^rf^ 



(1) 



where Hq is the Hubble constant, Qm is the density of mat- 
ter, c is the speed of light, L is the length of the box, and 
X are comoving distances where xo is the distance to the 
source galaxies. The summation is performed over the i* 
box. The number of particles associated with a pixel of the 
simulation is rip, the total number of particles within a sim- 
ulation is Nt, and s = Lp/L, where Lp is the length of the 
plane responsible for the lensing. R is the size of the 2D maps 
and Ar/. = rs^n^ where ri and r2 are comoving distances. 
For each of our 5 models, we have run 21 N-body simula- 
tions, each containing 256^ particles. We have used these 3D 
N-body simulations to derive 100 realizations of the conver- 
gence field for each model. The field is 3.95°x 3.95° and it has 
been downsampled to 1024 x 1024 pixels (1 pixel = 0.23^). 
Figure [2] shows one of these convergence maps. As said pre- 
viously, the 5 cosmological models have been selected along 
the ag-^m degeneracy corresponding to current constraints 
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from a power spectrum analysis. Additionally, we have ver- 
ified that the power spectrum of the 5 cosmological models 
are still degenerated with the survey area considered. 

The convergence map n that corresponds to the pro- 
jected matter density is not directly observable but, it can 
be derived from the observed shear maps 7 using the follow- 
ing relation ( Kaiser and Squires,1993^ ,Starck et al.||2006 ): 

^ = Pi7i + P272, (2) 



0.0 0.2 0.4 0.6 



with : 



A(k) 



kl 



2kik2 



A;2 



where the hat symbol denotes Fourier transform. Inversely, 
the shear maps 7^ can easily be derived from the convergence 
map using the following relation: 



7i = Pik 



(3) 



In practice, observed shear maps 7^^^* are obtained 
by averaging over a finite number of galaxies and are 
therefore noisy: 7°^* = 'ji -\- Ni, where A^i and N2 are 
noise contributions with zero mean and standard de- 
viation an = Gel \/ A.Ug with A being the area of the 
pixel in arcmin^ and Ug the average number of galaxies 
per arcmin^. Typical values for the density of galaxies 
is Ug — 30 gal/arcmin^ for ground-based simulations 
and Tig — 100 gal/arcmin^ for space-based simulations. 
Although a bit optimistic, these two configurations have 
been considered to derive noisy shear maps and to compute 
the shear statistics described in §3.1[ The derived noisy 
shear maps are downsampled to 1024 x 1024 pixels (1 
pixel = 0.23^) like the simulated convergence maps. An 
estimation of the corresponding noisy convergence maps 
can be derived from the equation 12]: kn — k -\- N, where 
N = PiNi + P2N2. As follows, the noise N in Kn is 
still Gaussian and uncorrelated. The inversion does not 
amplify the noise, but Kn may be dominated by the noise 
if N is large, which is the case in practice. Ground-based 
and space-based simulations of convergence maps have 
been derived by this way. The noisy convergence maps 
derived by inversion are still 3.95°x 3.95° downsampled to 
1024 X 1024 pixels (1 pixel = 0.23') and they have been 
used to compute the convergence statistics described in §3.3| 



3 WEAK LENSING STATISTICS 

The properties of the shear field 7 (or associated convergence 
field K,) can be measured statistically and reveal precious 
information about cosmological parameters. Up to now, 
cosmic shear studies have focused mainly on second-order 
statistics which only probe the Gaussian part of the matter 
distribution. However, the matter distribution is composed 
of no n- Gaussian features such as the clusters of galaxies 
that are the result of the non-linear evolution of the primor- 
dial Gaussian field. Therefore, higher-order statistics are re- 
quired to probe the non-Gaussian part of the field and thus 
improve our constraints on cosmological parameters. This 
study will focus on these higher-order statistics estimated 
both in the shear maps 7 and in the convergence maps k. 
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Figure 2. Simulated convergence map corresponding to a real- 
ization of a cosmological model with parameters: Qm = 0.23, 
n^ = 0.77, h = 0.594, erg = 1. The field is 3.95° x 3.95° down- 
sampled to 1024 X 1024 (pixel scale = 0.230- 



3.1 Shear statistics 

The shear field 7 can be directly derived from measure- 
ments of the shape of galaxies. For this reason, two-point 
statistics of the shear field have become the standard way 
of constraining cosmological parameters (see for example 
Maoh et aL][200T] [Hoekstra et ah] [2006) [Benjamin eFldl 
2007| |Fu et al.||2008t . Most of the interest in this type of 
analysis comes from its potential to constrain the spectrum 
of density fluctuations present in the late Universe and thus 
the cosmological parameters. However, as said previously, 
cosmological parameters cannot be determined accurately 
using only second-order statistics because only the Gaussian 
features of the field are captured by this method. Therefore, 
higher-order statistics have been introduced to probe the 
non- Gaussian features of the field and thus break degenera- 
cies. However, although the two-point correlations of the 
spin- 2 shear field 7^ can be reduced to a scalar quantity 
for parity reasons, this is not the case for higher-order 
moments of the shear field (see [Schneider and Lombardi 



2003 [Takada and Jain1[2003[ [Zaldarriaga and Scoccimarro 
2003). A way to get round this problem is to estimate 
higher-order statistics of the aperture mass Map, which has 



been introduced by Schneider (1996), rather than using the 
shear field directly. 

The aperture mass Map is one of the most widely 
used techniques for probing non-Gaussianity from the shear 
field (e.g. [Schneider et al.|1998[[Jarvis et al.|2004[[Kilbinger 



and Schneider[[2005[ [Dietrich and Hartlap[[2bl0 ) 

The aperture mass Map can be expressed in terms of 
the tangential component of the shear 7^ : 



Mapie) 



I' 



^'dQe{^)lt('d), 



(4) 



where Qe{^) is a radially symmetric, finite and continuous 
weight function and i9 is measured from the center of the 
aperture. The choice of the weight function Qe{'d) is arbi- 
trary at this point. In this paper, we have used the form 
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introduced by Crittenden et al. 2002, that has been found 
to be more sensitive for constraining Qm than other forms 
(Zhang et al, 2003) : 

"26'2 



9W 



^' 



47tO^ 



exp 



(5) 



In contrast with the shear field, which is a spin- 2 field from 
which higher-order moments are not trivial to define, the 
aperture mass defined by equation [4] is a scalar (spin-0) 
field, the skewness and kurtosis of which are well defined. 

In this section, for each noisy shear maps 7 described 
in q2] we have estimated the aperture mass Map{0) for 
several apertures from the relations p] and pj. Figure 
[3] shows some of these aperture mass maps for apertures 
e = 0.46', 0.92^ 1.85', 3.70', 7.40', 11.20'. Then, the following 
non- Gaussian statistics have been estimated: 

(i) The skewness of the aperture mass map (M^p) 

is the third-order moment of the aperture mass Map{0) and 
can be computed directly from shear maps filtered with 
different aperture mass. The skewness is a measure of the 
asymmetry of the probability distribution function. The 
probability distribution function will be more or less skewed 
positively depending on the abundance of dark matter 
haloes at the scale. The formalism exists to predict the 
skewness of the aperture mass map for a given cosmological 
model (see e.g. Jarvis et al.||2004" Kilbinger and Schneider 
[2005] ). 

(ii) The kurtosis of the aperture mass map {M^p) 

is the fourth-order moment of the aperture mass Map{0) 
and can be computed directly from the different aperture 
mass maps. The kurtosis is a measure of the peakedness of 
the probability distribution function. The presence of dark 
matter haloes at a given scale will flatten the probability 
distribution function and widen its shoulders leading to 
a larger kurtosis. The formalism exists to predict the 
kurtosis of the aperture mass map for a given cosmologi- 
cal model (Jarvis et al.|2004' Kilbinger and Schneider|2005 ) . 



(iii) The peak count of the aperture mass maps 

Pm^ ■ A peak is defined as connected pixels above a detec- 
tion threshold T. We consider all pixels that are connected 
via the sides or the corners of the pixel as one structure. 
It means, we are not discriminating between peaks due to 
massive halos and peaks due to projections of large-scale 
structures. The formalism exists to predict the peak counts 
in weak-lensing surveys, including the fraction of spurious 



detections caused by projections effects (e.g. Maturi et al, 
20T0I. 



We have reviewed the state-of-the-art of the non- 
Gaussian statistics used to constrain cosmology. Interesting 
analytical results relative to the shear three-point correla- 
tion function or the convergence bispectrum were also re- 
ported (e.g. Ma and Fry (2000a|b); Scoccimarro and Couch- 



(2001 ); Cooray and Hu (2001 )). However, when consid- 



ering only the equilateral configuration of the bispectrum, 
it has been shown that the discrimination efficiency of the 
cosmological models is relatively poor ( Pires et al.||2009b|). 
An analytical comparison has been performed in |Berge et al.] 
([2010| (for an Euclid-like survey) between the full bispec- 



trum and an optimal match-filter peak count and both ap- 
proaches were found to provide similar results. However, as 
the full bispectrum calculation has a much higher complex- 
ity than peak counting, and no public code exists to compute 
it (only an equilateral code is available (Pires et al. 2009a)), 
the bispectrum has not been considered in this study. 



3.2 Convergence statistics 

In this section, we have used the noisy convergence maps 
described in ^ The convergence has already been used in 
some studie s (^ernardeau et al.|1997b[ Hamana et al.|2004 



Pires et al.||2009a| |Wang et al.||2009| [Berge et al.||20iot To 
extract non-Gaussian information from higher-order statis- 
tics. In this paper, we want to study the ability of higher- 
order shear statistics compared to higher-order convergence 
statistics to break the ag-^m degeneracy. However, a fair 
comparison requires to compare the previous shear statis- 
tics with the convergence statistics at the same scale of 
the aperture mass. We could have used the definition of the 
aperture mass expressed in terms of the convergence given 
by: 



Map{e)= f d^i)UeWK{i)), 



with: 



f/eW = 



^' 



27r6'2 \ 26»2 



exp 



26*2 



(6) 



(7) 



However, we have preferred to use an undecimated isotropic 
wavelet transform: the "a trous" wavelet transform, which 
computes simultaneously J aperture mass maps for dyadic 
scales. The "a trous" wavelet transform decomposes a 

convergence map k (of size A^ x A^) as a superposition of the 
form: 






(8) 



The algorithm outputs J + 1 sub-band arrays of size N x N 
where cj is a coarse or smooth version of the original ima ge k 
and Wj represents the details of k at scale 2^ ( see Starck et al. 
1998; Starck and Murtagh||2002| for details) . [Leonard et al 



(2011) have shown that the wavelet bands Wj are formally 
identical to aperture mass maps at scale — 2^ except the 
filter Ue is replaced by the following wavelet function ip{x,y): 






(9) 



where (p{x,y) = (p(x)(p{y) and (p{x) is a compact function 
(a B3-spline function) defined by: 



(p{x) 



12 



(|x-2|^-4|x- 



11' + 61x1^ 



-4|x + l|' + |x + 2|^) 



Fig. |4] displays in the Fourier domain in solid lines the aper- 
ture mass filters at scale ^^ = 2, 4, 8, 16, 32 pixels (1 pixel = 
0.23 arcmin) and in dashed lines, the corresponding wavelet 
filters at the same scale 2^ pixels with j — 1, 2, 3, 4, 5. To as- 
sess the response of aperture mass filters and wavelet filters, 
we have generated artificial shear data from a null conver- 
gence map with a single central delta function. The aperture 
mass algorithm has then been applied to the resulting shear 
maps and the wavelet transform has been computed from 
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Figure 3. Aperture mass maps obtained from noisy shear maps (space-based simulations) filtered with an aperture mass with scales 
= 0.46^ 0.92^ 1.85^ 3.70^ 7.40^ 11.20^ The field is 3.95° x 3.95° downsampled to 1024 x 1024 i.e. a pixel scale of 0.23^ 
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Figure 4. Frequency response of the aperture mass filters for 
scales ei= 0.46^ 0.92^ 1.85^ 3.70^ 7.40' (solid lines) and fre- 
quency response of the wavelet filters at the same scales (dashed 
lines). 



the convergence map. Then the response of these filters in 
Fourier space has been obtained by considering their power 
spectra. 

These two filter banks are really close, except the last 
filters (from the left). The last wavelet filter is a high-pass fil- 
ter whereas the last aperture mass filter is a band-pass filter. 
As a consequence, we can access to the finest scales of the 
image using the wavelet transform method which is not the 
case with the aperture mass method. Nevertheless, the gen- 
eralized definition of the aperture mass statistics estimated 



from the shear n-point correlation functions (see Schneider 
et al.|2005 Semboloni et al.|2011) also gives access to these 



small scales because it is estimated on the shear catalogue 
directly. However, this has not been considered in this study 
because it is too intensive computationally. In the same way, 
there exists also some wavelet transforms that can work di- 
rectly on the shear catalogue (see e.g. Deriaz et al., 2012). 
This has not been considered in this analysis because the 
resolution of the maps is very good and the noise is already 
dominant at this scale. Considering now the other filters, 
we see that the aperture mass filters have some unwanted 
oscillations that make the wavelet filters much more local- 
ized in Fourier space. These oscillations in Fourier space are 
due to the fact that the aperture mass Qei^i^) defined in [s] 
is truncated for i^ > (see [Leonard et al.|[2011| for more 
details). Another point in favour of the wavelet transform 
is its time to compute. The wavelet transform complexity is 
oc 0{N'^{j + 1)) compared to the aperture mass complexity, 
which is oc 0{N'^ ^- Of). For one of our weak lensing simu- 
lations (of size 1024 x 1024), it is about 250 times faster to 
compute the 5 considered scales with a wavelet transform 
than with the aperture mass definition. The comparison be- 
tween the aperture mass and the wavelet transform method 
has only been performed for scales : 0.92', 1.85', 3.70' and 
7.40'. The scales in between have not been considered in 
this analysis because the dyadic scales are sufficient to char- 
acterize the variation of the discrimination as a function of 
the scale. Nevertheless, if we want to perform a more precise 



analysis, the scales in between this dyadic scaling can easily 
be obtained with the wavelet transform method by chang- 
ing the initial pixel scale of the map. In this study, for each 
simulated noisy convergence map n described in ^ we have 
estimated the wavelet transform from the previous definition 
[s] and then the following statistics have been estimated for 
each wavelet band Wj: 

(i) The skewness of the wavelet band (w^) that is com- 
puted directly from the different wavelet bands. 

(ii) The kurtosis of the wavelet band (wj) that is com- 
puted directly from the different wavelet bands. 

(iii) The peak count of the wavelet band P^. . A peak is 
defined as connected pixels above a detection threshold T. 



3.3 Convergence statistics in denoised maps 

In this section, we want to derive higher-order statistics 
from denoised maps because the more the data are noisy 
and the more the probability distribution function looks like 
a Gaussian, the less the higher-order statistics will be useful. 
We expect, for example, the skewness and the kurtosis 
to tend to zero with an additive Gaussian noise and the 
clusters to be more difficult to extract. Therefore, to extract 
the non-Gaussian structures and reduce the impact of the 
noise in the analysis, we have used the MRLens denoising 
proposed by Starck et al. (2006), which is a multiscale 



Bayesian denoising based on the sparse representation of 
the data. In [Starck et aLJ (|2006t; [Teyssier et~^ (|2009|); 



Pires et al. (2009a), the authors have shown that this 



method outperforms several standard techniques to detect 
non-Gaussian structures such as Gaussian filtering, Wiener 
filtering and MEM filtering. 

In this study, the MRLens denoising has been used to 
denoise each of the simulated convergence maps k. The MR- 
Lens denoising software is available at the following address: 
"http://irfu.cea.fr/Ast/mrlens_software.php". We have only 
applied this denoising to the noisy convergence maps be- 
cause of the difficultly of applying denoising to the spin-2 
shear field. Then, a wavelet transform has been applied to 
the denoised convergence map k in order to compute the de- 
noised wavelet bands Wj and estimate the following statis- 
tics: 

(i) The skewness of the denoised wavelet band (w^) that 
is computed directly from the denoised wavelet band j. 

(ii) The kurtosis of the denoised wavelet band {wj) that 
is computed directly from the denoised wavelet band j. 

(iii) The peak count of the denoised wavelet band Pj,. . 
A peak is defined as connected pixels above a threshold e. 
Contrary to §3.1| and §3.2| where the threshold T is used to 
extract peaks from the noise, the threshold e (set to the rms 
value of the denoised wavelet band Wj) is only used to reject 
spurious detections in the denoised convergence maps. 



4 RESULTS 

4.1 The methodology 

In this study, we are interested in comparing the ability of 
the previous statistics to break the ag-^m degeneracy, using 
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92.35 % 


85.45 % 








1.85' 


03.45 % 


01.45 % 


91.25 % 


89.20 % 


2.80' 


03.40 % 


07.95 % 


83.25 % 


90.00 % 


3.70' 


15.15 % 


23.00 % 


69.40 % 


86.70 % 


5.60' 


28.55 % 


28.90 % 


02.45 % 


70.55 % 


7.40' 


26.95 % 


24.30 % 


4.90 % 


60.50 % 


11.20' 


23.00 % 


20.30 % 


09.65 % 


01.90 % 



^. 


{Mapio.r) 


(M.^(^,)4> 


p2cr 


p3cr 


0.46' 


01.95 % 


02.25 % 


03.45 % 


06.50 % 


0.69' 


01.75 % 


01.05 % 


12.50 % 


05.60 % 


0.92' 


02.65 % 


01.30 % 


31.95 % 


21.90 % 


1.40' 


07.70 % 


01.75 % 


49.35 % 


41.90 % 


1.85' 


04.05 % 


02.15 % 


55.85 % 


48.85 % 




01.80 % 


03.15 % 


55.00 % 






2.80' 


62.65 % 










3.70' 


04.00 % 


05.20 % 


54.40 % 


54.45 % 


5.60' 


09.70 % 


10.90 % 


11.05 % 


60.40 % 


7.40' 


12.00 % 


10.25 % 


14.15 % 


51.90 % 


11.20' 


09.05 % 


08.90 % 


04.30 % 


09.60 % 



Table 1. Mean discrimination efficiency (in percent) from noisy 
aperture mass maps (space-based simulations) for scales 6i = 
0.46^ 0.92^ 1.40^ 1.85^ 2.80^ 3.70^ 5.60^ 7.40^ 11.20^ 



Table 2. Mean discrimination efficiency (in percent) from noisy 
aperture mass maps (ground-based simulations) for scales 6i = 
0.46^ 0.92^ 1.40^ 1.85^ 2.80^ 3.70^ 5.60^ 7.40^ 11.20^ 



the set of cosmological simulations described in ^ For this 
purpose, we have characterized, for each statistic, its abilty 
to discriminate between the 5 different cosmological models. 

(f2009a), we have computed 



As has been done in Pires et al. 



a "discrimination efficiency" that expresses in percentage 
the ability of a statistic to discriminate between two cos- 
mological models. For this purpose, a statistical tool called 
FDR (False Discovery Rate) introduced by Benjamin! and 
Hochberg ( 1995) has been used to set in an adaptive way 
the thresholds to classify between the different cosmologi- 
cal models. Each threshold is estimated in such way that 
the rate of allowed false detections is inferior to a 0.05. The 
larger the discrimination efficiency is, the less the proba- 
bility distributions of the statistic values for the different 
cosmological models overlap. The optimal statistic will be 
the one that maximizes the discrimination for all pairs of 
models. A mean discrimination efficiency can be estimated 
for each statistic by averaging the discrimination efficiency 
across all the pairs of models. 



4.2 Shear statistics results 

Tableland Table [2] show the mean discrimination efficiency 
for the shear statistics, described in p.l| estimated for vari- 
ous aperture mass Map{Oi) respectively for space-based and 
ground-based simulations. 

On first glance at these two tables, we can see that 
the results worsen with noise whatever the statistics. As ex- 
pected, the skewness and the kurtosis are very poor at dis- 
criminating between different cosmological models in noisy 
aperture mass maps because the aperture mass map prob- 
ability distribution function tends to a Gaussian distribu- 
tion as the noise increases. This makes the skewness and the 
kurtosis of the aperture mass map tend to zero. The result 
with peak counting is significantly better because a basic de- 
noising is applied by only selecting the peaks above a given 
threshold in the aperture mass map. However, the choice of 
this threshold is relatively important because the differences 
between a 2a and a 3a threshold are significant. 





Spc 


ice-based 


Ground-based 




Purity 


Completeness 


Purity 


Completeness 


p2cr 

^Map 


16.07 % 


61.27 % 


10.77 % 


35.60 % 


p3cr 

^Map 


47.63 % 


35.11 % 


33.15 % 


11.56 % 



Table 3. Purity and completeness for the peak count for space- 
based (left) and ground based (right) aperture mass maps corre- 
sponding to realizations of the cosmological model with Qm = 0.3 
and (78 = 0.9). P^ (respectively P^ ) is defined for peaks 
above a 2cr-threshold (respectively a 3cr-threshold) on noisy aper- 
ture mass maps for scale = 1.85^ 



Table [3] shows the purity and the completeness for 
the peak count estimated on noisy aperture mass maps 
{0 = 1.85^) for space-based (left) and ground based (right). 
Completeness and purity are two important criteria to eval- 
uate the performance of a peak detection method. Purity is 
defined as the ratio of true detections to the total number 
of peaks detected, and completeness is defined as the ratio 
of true detections to the total number of peaks in the sim- 
ulation. The total number of peaks per scale is estimated 
from the simulated 2D convergence maps (without noise) 
with different aperture masses. It means, we are not discrim- 
inating between peaks due to massive halos and peaks due 
to projections of large scale-structures. With a 2cr-threshold, 
the completeness is maximal but the purity is poor because 
there is a large number of false detections due to shot noise, 
among the total number of detected peaks. Thus, the num- 
ber of detected peaks is considerably overestimated espe- 
cially at small scales for which the noise is important. Thus, 
the choice of the best threshold is a trade-off between purity 
and completeness. 

Another important parameter is the scale 6 of the aper- 
ture mass. The discrimination efficiency depends strongly 
on the scale that is considered. The best discrimination effi- 
ciency scale is displayed in bold, for each statistic, in Table 
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Scale 


{wp 


{wp 


p2cr 


p3cr 


Finest scales 


02.00 % 


01.15 % 


12.05% 


00.70 % 


0.92' 


37.95 % 


04.75 % 


86.30 % 


73.05 % 




03.55 % 


02.10 % 






1.85' 


94.40 % 


93.85 % 








3.70' 


18.25 % 


25.65 % 


84.05 % 


87.05 % 


7.40' 


36.40 % 


30.90 % 


24.60 % 


66.35 % 



Table 4. Mean discrimination efficiency (in percent) from noisy 
convergence wavelet maps (space-based simulations) for the finest 
scales of the image given by the high-pass filter (first column) and 
for scales 0.92^ 1.85^ 3.70^ 7.40' (other columns). 



Table 5. Mean discrimination efficiency (in percent) from noisy 
convergence wavelet maps (ground-based simulations) for the 
finest scales of the image given by the high-pass filter (first col- 
umn) and for scales 0.92^ 1.85^ 3.70^ 7.40' (other columns). 



[T] and Table [2] We see that the best scale depends on the 
statistic that is used. This difference can be explained by 
the fact that the statistics are not sensitive in the same way 
to the different characteristics of the clusters. Skewness and 
kurtosis are very sensitive to the density of the clusters, 
e.g. very dense clusters will skew significantly the probabil- 
ity distribution function whereas small clusters will have a 
small impact. In contrast, the peak count is mainly sensitive 
to the number of clusters regardless of their masses or their 
density. A massive cluster will be accounted in the same way 
as a small cluster if it is detected. 

The best discrimination efficiency in noisy aperture 
mass maps (for space-based simulations) has been obtained 
with the peak count with a 2a threshold (92.35 %) for a 
scale of 1.40^ This is definitely better than the best result 
obtained with the skewness (33.40 %) as well as the best 
result obtained with the kurtosis (28.90 %). 



4.3 Convergence statistics results 

Table [4] and Table [5] show the mean discrimination efficiency 
for the convergence statistics, described in §3.2| estimated 
for different wavelet scales (2-^) respectively for space-based 
and ground-based simulations. As previously, the skewness 
and the kurtosis are very poor at discriminating between 
different cosmological models in noisy convergence maps be- 
cause the skewness and the kurtosis tend to zero as noise is 
increased (see Fig.lsl). 

A comparison with Table IT] and Table [2] shows that the 





Space-based 


Ground-based 




Purity 


Completeness 


Purity 


Completeness 


p2cr 


22.24 % 


63.20 % 


14.85 % 


36.31 % 


p3cr 

^W3 


56.66 % 


38.76 % 


42.02 % 


13.11 % 


pmrlens 


84.30 % 


49.55 % 


75.37 % 


25.92 % 



Table 6. Purity and completeness for the peak count for space- 
based (left) and ground-based (right) convergence maps corre- 
sponding to realizations of the cosmological model with Qm = 0.3 
and as = 0.9. P^^ (respectively P^^) is defined for peaks above a 
2cr-threshold (respectively a 3cr-threshold) on noisy convergence 
maps at the third scale of a wavelet transform (1.85') and p^^^^^^ 
is defined for peaks above a e-threshold on MRLens denoised con- 
vergence maps at the third scale of a wavelet transform (1.85'). 



Scale 


(^|> 


(%"> 


p2cr 


p3cr 


Finest scales 


00.60 % 


00.65 % 


01.20 % 


00.65 % 


0.92' 


07.50 % 


01.60 % 


38.15 % 


13.45 % 




06.70 % 


02.25 % 






1.85' 


62.75 % 


58.85 % 








3.70' 


04.15 % 


06.45 % 


62.25 % 


52.40 % 


7.40' 


14.55 % 


17.65 % 


44.20 % 


40.90 % 



results obtained with the aperture mass maps at the same 
scale are very similar. This is not a surprising result be- 
cause it has been shown by [Leonard et al." (2011) that ap- 
plying aperture mass filters at dyadic scales in shear maps 
is comparable to performing a wavelet transform of the con- 
vergence map that can be derived from the shear maps by 
inversion in Fourier space. This also explains the similarity 
of the wavelet filters compared to the aperture mass filters at 
the same scales in Fig. [4] However, the results are slightly im- 
proved with the wavelet transform for every statistic and ev- 
ery scale, which tends to show that the shape of the wavelet 
filters is more efficient to capture the non-Gaussian struc- 
tures present in the weak lensing maps. This is possibly a 
consequence of the oscillations seen in the aperture mass fil- 
ters (see Fig. p|, which gives rise to a small leakage of the 
signal into higher frequencies. 

Some other studies have been conducted to find 
an optimal filter for detecting dark matter haloes 
dMaturi et al.|2QQ5'||Pace et a l. 2007) and thus avoid 
the spurious peaks due to large-scale structure pro- 
jections. However, these filters are less efficient be- 
cause the projection effects that are normally a main 
source of uncertainty Avhen probing the clusters, 
here serve as an additional source of cosmological 
information (see Dietrich and Hartlap 2010[ Wang 
et al.||2009| for more details) 



The best discrimination efficiency in noisy convergence 
maps (for space-based simulations) has been obtained with 
the peak count (94.40 %), for a scale of 1.85^ and a 2cr- 
threshold. 

Table [6] shows the purity and the completeness for the 
peak count at the third scale of a wavelet transform (1.85') 
for space-based (left) and ground-based (right) convergence 
maps. A comparison with Table [3] shows that both purity 
and completeness are improved with the wavelet transform. 
As with the aperture mass statistic, the completeness is 
maximal with a 2cr-threshold. 

As previously, the constraints on cosmological models 
obtained with peak count (94.40 %) are significantly better 
than the ones that can be reached with the skewness (37.95 
%) and the kurtosis (30.90 %). 
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Scale 


(^f> 


{wp 


P^, 


Finest scales 


53.40 % 


43.20 % 


68.35 % 


0.92' 


47.90 % 


41.15 % 


92.45 % 




58.80 % 


44.70 % 






1.85' 


96.75 % 










3.70' 


63.30 % 


48.05 % 


90.40 % 


7.40' 


54.90 % 


40.45 % 


63.45 % 



Scale (j) 
4 



Table 7. Mean discrimination efficiency (in percent) from MR- 
Lens denoised convergence wavelet maps (space-based simula- 
tions) for the finest scales of the image given by the high-pass 
filter (first column) and for scales 0.92^ 1.85^ 3.70^ 7.40^ (other 
columns). 



Scale 


(^|> 


{wp 


P^, 


Finest scales 


42.15 % 


30.05 % 


38.20 % 


0.92' 


35.95 % 


28.60 % 


40.45 % 


1.85' 


31.65 % 


20.85 % 


62.35 % 




41.80 % 


29.95 % 






3.70' 


72.65 % 










7.40' 


44.75 % 


32.25 % 


54.55 % 



Table 8. Mean discrimination efficiency (in percent) from MR- 
Lens denoised convergence wavelet maps (ground-based simula- 
tions) for the finest scales of the image given by the high-pass 
filter (first column) and for scales 0.92^ 1.85^ 3.70^ 7.40^ (other 
columns). 



4.4 Denoised Convergence statistics results 

In this section, we want to show that the convergence statis- 
tics can be improved significantly if denoising is applied to 
the convergence maps. As said previously, the convergence 
maps have been denoised using the MRLens denoising de- 
scribed in |Starcket"aL| ( |2006t . Table |7| and Table [s] show the 
mean discrimination efficiency for the denoised convergence 
statistics, described in §3.3[ estimated for different wavelet 
scales (2-^) respectively for space-based and ground-based 
simulations. 

As expected, the MRLens denoising improves consider- 
ably the discrimination efficiency of the skewness and kur- 
tosis. This comes from its ability to reconstruct the non- 
Gaussian structures that dominate at small scales. However, 
the skewness and kurtosis values are significantly overesti- 
mated compared to original kurtosis, as shown in Fig. |5] 
because the MRLens denoising is only efficient in recover- 
ing high peaks in the signal, which affects the tails of the 
probability distribution function. 

The MRLens denoising also improves the discrimina- 
tion efficiency of the peak count at all scales, especially for 
ground-based simulations for which the noise is important. 

Table [6] shows the purity and the completeness for 
peak counting on MRLens denoised maps for space-based 
and ground-based simulations. In the MRLens denoising, we 
again have the usual trade-off between purity and complete- 
ness. A different threshold is selected for each wavelet band. 
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Figure 5. Top: Mean skewness per scale for original convergence 
maps (black) , space-based noisy convergence maps (red) and MR- 
Lens denoised convergence maps (blue). Bottom: Mean kurtosis 
per scale for original convergence maps (black) , space-based noisy 
convergence maps (red) and MRLens denoised convergence maps 
(blue). The skewness and the kurtosis of the noisy convergence 
maps are considerably reduced especially at small scales for which 
the noise is important. In contrast, the skewness and the kurtosis 
are significantly overestimated on MRLens denoised convergence 
maps. These convergence maps correspond to realizations of the 
cosmological model with Qrn = 0.3 and (jg = 0.9). 



and this is done in an adaptive way, conformed to a False 
Discovery Rate method ( |Starck et al.|2006 ), which provides 
a more robust discrimination. The completeness with MR- 
Lens denoising is slightly inferior to a 2cr-threshold but its 
purity is maximal. 

The best discrimination efficiency in denoised conver- 
gence maps (for space-based simulations) has been obtained 
with the peak count (96.75 %) still for a scale of 1.85^, in 
perfect agreement with the results of Fires et al. (2009a). 



Table [9] shows the discrimination efficiency obtained with 
this statistic that enables to discriminate between the five 
cosmological models even for contiguous models for which 
the discrimination is challenging. The Table is not symmet- 
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model 1 


model 2 


model 3 


model 4 


model 5 


model 1 


X 


85 % 


100 % 


100 % 


100 % 


model 2 


89% 


X 


92 % 


100 % 


100 % 


model 3 


100 % 


92 % 


X 


89 % 


100 % 


model 4 


100% 


100% 


92% 


X 


98% 


model 5 


100 % 


100 % 


100 % 


98 % 


X 



Table 9. Discrimination efficiency (in percent) between the 5 cosmological models obtained with the peak count on denoised convergence 
maps at the third scale (1.85') of a wavelet transform (space-based simulations). 



ric because the probability distributions of the statistics are 
not symmetric and to quantify the discrimination between 
two cosmological models, the FDR method sets two different 
thresholds in an adaptive way. This result can be compared 
with the result obtained in Pires et al. (2009a') (Table 7) with 
another set of cosmological simulations. The results are very 
similar. The small discrepancies are only due to the limited 
size of the cosmological simulation sample. 



5 CONCLUSION 

The goal of this paper was to investigate how to best extract 
non-Gaussianity from weak lensing surveys to constrain the 
cosmological model. For this purpose, we have been inter- 
ested in showing that there is an extra information that can 
be derived from higher-order statistical analysis of the con- 
vergence maps beyond what can be obtained from higher- 
order statistical analysis of the shear maps. Therefore, we 
have compared the efficiency of several higher-order shear 
and convergence statistics to break the ag — ^m degener- 
acy, by comparing their ability to discriminate between 5 
cosmological models along this degeneracy. 

Most of the techniques used to estimate higher-order 
statistics from the spin-2 shear field are based on the aper- 
ture mass expressed in terms of the tangential component of 
the shear (see relation [4]). Analogous convergence statistics 
can be obtained by using the aperture mass defined from 
the convergence maps (see relation [6]). However, in accor- 
dance with Leonard et al.| ( |2011| , we have preferred to use 
an alternative solution that computes simultaneously mul- 
tiple aperture mass maps for dyadic scales: the "a trous" 
wavelet transform. In this study, we have observed that this 
method is 250 times faster than the aperture mass method 
to run on our simulations and that the wavelet filters are 
much more localized in Fourier space compared to aperture 
mass filters. It follows that the results obtained in wavelet 
convergence maps compared to the results in aperture mass 
maps are very similar but slightly improved in the wavelet 
case for every statistic and every scale. Therefore, contrary 
to a generally accepted idea, the noise properties in aperture 
mass maps are not better than in convergence maps, if the 
same scale is considered. 

Contrary to another accepted idea, further important 
cosmological information can be extracted from noisy con- 
vergence maps if a denoising such as MRLens is used. This 
comes from its ability to reconstruct the non-Gaussian struc- 
tures that are induced by the growth of structures. In this 



study, we have shown that the MRLens denoising improves 
considerably the discrimination efficiency of the skewness 
and kurtosis. It also improves the discrimination efficiency 
of the peak count especially for ground-based simulations 
(n^ = 30 gal/arcmin^) for which the noise is important. 

For an Euclid-like survey, the density of galaxies is ex- 
pected to be around Ug = 40 gal/arcmin^ for the wide-field 
survey and around Ug — 80 gal/arcmin^ for the deep-field 
survey. It is clear from this study, that the non-Gaussian 
statistical analysis should be performed in denoised conver- 
gence maps as described in §3.3| 

Finally, the best non-Gaussian statistic to constrain cos- 
mological model in combination with the power spectrum 
has been found to be the peak count per scale. And further 
cosmological information should be obtained by combining 
the constraints obtained with the peak count at different 



scale as shown by Marian et al. (2011 ). This will be investi- 



gated in a future work. 



ACKNOWLEDGMENTS 



This work has been supported by the European Research 
Council grant SparseAstro (ERC-228261). 



REFERENCES 

Bacon, D., Refregier, A., and Ellis, R.: 2000, MNRAS ^^ 

318-625 
Bacon, D. J., Massey, R. J., Refregier, A. R., and Ellis, 

R. S.: 2003, MNRAS3AA, 673 
Bartelmann, M. and Schneider, P.: 2001, Phys. Rep. 340, 

291 
Benjamin, J., Heymans, C, Semboloni, E., van Waerbeke, 

L., Hoekstra, H., Erben, T., Gladders, M. D., Hetter- 

scheidt, M., Mellier, Y., and Yee, H. K. C: 2007, MNRAS 

381, 702 
Benjamini, Y. and Hochberg, Y.: 1995, J. R. Stat. Soc. B 

57, 289 
Berge, J., Amara, A., and Refregier, A.: 2010, ApJ 712, 

992 
Bernardeau, F., van Waerbeke, L., and Mellier, Y.: 1997a, 

A&A 322, 1 
Bernardeau, F., van Waerbeke, L., and Mellier, Y.: 1997b, 

A&A 322, 1 
Cooray, A. and Hu, W.: 2001, Aj9j 548, 7 
Dahle, H.: 2006, A^^J 653, 954 



Cosmological constraints from the capture of non-Gaussianity in Weak Lensing data 11 



Ellis, R., and 



and Bartelmann, M.: 



Dietrich, J. P. and Hartlap, J.: 2010, MNRAS 402, 1049 
Fu, L., Semboloni, E., Hoekstra, H., Kilbinger, M., van 

Waerbeke, L., Tereno, L, Mellier, Y., Heymans, C, 

Coupon, J., Benabed, K., Benjamin, J., Bertin, E., Dore, 

O., Hudson, M. J., Ilbert, O., Maoli, R., Marmo, C, Mc- 

Cracken, H. J., and Menard, B.: 2008, A&A 479, 9 
Hamana, T., Takada, M., and Yoshida, N.: 2004, MNRAS 

350, 893 
Hoekstra, H. and Jain, B.: 2008, Annual Review of Nuclear 

and Particle Science 58, 99 
Hoekstra, H., Mellier, Y., van Waerbeke, L., Semboloni, 

E., Fu, L., Hudson, M. J., Parker, L. C, Tereno, L, and 

Benabed, K.: 2006, ApJ647, 116 
Jarvis, M., Bernstein, G., and Jain, B.: 2004, MNRAS 362, 

338 
Kaiser, N. and Squires, G.: 1993, Ap J 404, 441 
Kaiser, N., Wilson, G., and Luppino, G.: 2000, AJpp 318- 

625 
Kilbinger, M. and Schneider, P.: 2005, A&A 442, 69 
Leonard, A., Pires, S., and Starck, J.-L.: 2011, MNRAS 
Ma, C.-P. and Fry, J. N.: 2000a, Ap J 643, 503 
Ma, C.-P. and Fry, J. N.: 2000b, Ap J 638, L107 
Maoh, R., Van Waerbeke, L., Mellier, Y., Schneider, P., 

Jain, B., Bernardeau, F., Erben, T., and Fort, B.: 2001, 

A&A 368, 766 
Marian, L., Smith, R. E., Hilbert, S., and Schneider, P.: 

2011, ArXiv e-prints 
Massey, R., Refregier, A., Bacon, D. J. 

Brown, M. L.: 2005, MNRAS 369, 1277 
Maturi, M., Angrick, C, Pace, F., 

2010, A&A 519, A23 
Maturi, M., Meneghetti, M., Bartelmann, M., Dolag, K., 

and Moscardini, L.: 2005, A&A 442, 851 
Mellier, Y.: 1999, ARA&A 37, 127 
Munshi, D., Valageas, P., van Waerbeke, L., and Heavens, 

A.: 2008, Phys. Rep. 462, 67 
Pace, F., Maturi, M., Meneghetti, M., Bartelmann, M., 

Moscardini, L., and Dolag, K.: 2007, A&A 471, 731 
Pires, S., Starck, J.-L., Amara, A., Refregier, A., and 

Teyssier, R.: 2009a, A&A 505, 969 
Pires, S., Starck, J.-L., Amara, A., Teyssier, R., Refregier, 

A., and FadiU, J.: 2009b, MNRAS 396, 1265 
Refregier, A.: 2003, ARA&A 41, 645 
Refregier, A., Rhodes, J., and Groth, E. J.: 2002, Ap J 672, 

L131 
Schneider, P.: 1996, MNRAS 283, 837 
Schneider, P., Kilbinger, M., and Lombardi, M.: 2005, A&A 

431, 9 
Schneider, P. and Lombardi, M.: 2003, A&A 397, 809 
Schneider, P., van Waerbeke, L., Jain, B., and Kruse, G.: 

1998, MNRAS 296, 873 
Scoccimarro, R. and Couchman, H. M. P.: 2001, MNRAS 

325, 1312 
Semboloni, E., Schrabback, T., van Waerbeke, L., Vafaei, 

S., Hartlap, J., and Hilbert, S.: 2011, MNRAS 410, 143 
Starck, J.-L. and Murtagh, F.: 2002, Astronomical Image 

and Data Analysis, Springer- Verlag 
Starck, J.-L., Murtagh, F., and Bijaoui, A.: 1998, Image 

Processing and Data Analysis: The Multiscale Approach, 

Cambridge University Press 
Starck, J.-L., Pires, S., and Refregier, A.: 2006, A&A 451, 

1139 



Takada, M. and Jain, B.: 2003, MNRAS 344, 857 
Takada, M. and Jain, B.: 2004, MNRAS 348, 897 
Teyssier, R.: 2002, A&A 385, 337 
Teyssier, R., Pires, S., Prunet, S., Aubert, D., Pichon, C, 

Amara, A., Benabed, K., Colombi, S., Refregier, A., and 

Starck, J.-L.: 2009, A&A 497, 335 
Vale, C. and White, M.: 2003, Ap J 692, 699 
Van Waerbeke, L., Mellier, Y., Erben, T., Cuillandre, J., 

and Bernardeau, F. e. a.: 2000, A&A p. 318:30 
Wang, S., Haiman, Z., and May, M.: 2009, Ap J 691, 547 
Wittman, D., Tyson, J., Kirkman, D.and DellAntonio, L, 

and Bernstein, G.: 2000, Nature p. 405:143 
Zaldarriaga, M. and Scoccimarro, R.: 2003, ApJ 684, 559 



